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The role of interchain hopping in quasi-one-dimensional (Q-ID) electron systems is investi- 
gated by extending the Kadanoflt- Wilson renormalization group of one-dimensional (ID) sys- 
tems to Q-ID systems. This scheme is applied to the extended Hubbard model to calculate the 
temperature (T) dependence of the magnetic susceptibility, x{T). The calculation is performed 
by taking into account not only the logarithmic Cooper and Peierls channels, but also the non- 
logarithmic Landau and finite momentum Cooper channels, which give relevant contributions 
to the uniform response at finite temperatures. It is shown that the interchain hopping, t±, 
reduces xiT) at low temperatures, while it enhances x{T) at high temperatures. This notable 
dependence is ascribed to the fact that t±^ enhances the antiferromagnetic spin fluctuation 
at low temperatures, while it suppresses the ID fluctuation at high temperatures. The result is 
at variance with the random-phase-approximation approach, which predicts an enhancement of 
X(T) by over the whole temperature range. The influence of both the long-range repulsion 
and the nesting deviations on x{T) is further investigated. We discuss the present results in 
connection with the data of xiT) in the (TMTTF)2X and (TMTSF)2X series of Q-ID or- 
ganic conductors, and propose a theoretical prediction for the effect of pressure on magnetic 
susceptibility. 

KEYWORDS: quasi-one dimension, magnetic susceptibility, renormalization group, extended-Hubbard 
model, long-range interaction, nesting deviation, crossover, (TMTTF)2X, (TMTSF)2X 



1. Introduction 

Interacting electrons in ordinary metals and semicon- 
ductors are usually described in the framework of the 
Landau's Fermi liquid (FL) theory.^ What lies at the root 
of the Fermi liquid theory is the concept of single-particle 
excitations that behave as effective non-interacting elec- 
trons at low-energy. In one dimension, however, the 
FL picture is no longer valid due to the strongly en- 
hanced correlation effects, which lead to the absence 
of quasi-particles and the emergence of decoupled col- 
lective spin and charge excitations, which characterize 
different quantum liquid phases such as the Tomonaga- 
Luttinger liquid (TLL).^ While properties of the FL and 
TLL are independently well characterized, their relative 
importance for weakly coupled chains in the quasi-one- 
dimensional (Q-ID) case is much less understood. In Q- 
ID systems, a dimensional crossover between TLL and 
FL phases is expected to take place on the temperature 
scale. How such a crossover is achieved in practice is 
an important issue in weakly coupled chains of organic 
conductors (TMTCF)2X (C=S, Se, FFg, Br, CIO4, 
AsFe, etc.). 

Such a crossover has been observed in transport prop- 
erties of these materials.^ The temperature (T) depen- 
dence of resistivity (p) exhibits a power-law dependence 
p ^ T" (a ~ 0.5) for T above - lOOK, which is 
consistent with a TLL behavior.^ It is followed by a 
dependence below Tx, which is expected in a FL. For 
the uniform magnetic susceptibility, Xi however, the ex- 
pected temperature dependence below is not known 



theoretically. On experimental grounds, x(T) does not 
show any change around Tx,**'^ and is still temperature 
dependent. This indicates that the influence of the di- 
mensional crossover may be quite specific to observable 
quantities. It is of particular interest to examine how the 
physical quantities of the Q-ID correlated electron sys- 
tems are modified through T^. 

Some basic aspects of this dimensional crossover can be 
easily illustrated by first looking at the non-interacting 
case. When the interchain hopping tx is much smaller 
than the intrachain one t, the Fermi surface consists of 
two slightly warped planes as shown in Fig. 1, where the 
energy scale of the warping is of order of t± . Now when 
the temperature scale is much higher than t± , the warp- 
ing is incoherent due to thermal fluctuations, the quan- 
tum mechanical coherence of electrons is then thermally 
confined along the chains and the system is effectively 
one dimensional (ID). By contrast, for the energy scale 
much lower than tj^, the quantum coherence of electron 
extends over several chains, the warping of the Fermi sur- 
face becomes effective, and the system is considered as 
two- or three-dimensional. The scale t± thus corresponds 
to the crossover temperature Tx for non-interacting elec- 
trons. The presence of interactions, however, modifies 
this picture by reducing the crossover scale, which leads 
T* - tit^/ty/'^ with < (/)< 1.^ 

The main theoretical difficulty in studying the Q-ID 
systems resides in the treatment of the ID fiuctuations 
which affect the coupling between chains and the restora- 
tion of a FL behavior. At the lowest order of ID per- 
turbation theory, non-FL behavior basically comes from 
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the quantum interference between the Cooper and the 
Peierls infrared instabiUties. This is an effect that is not 
captured by approaches like random phase approxima- 
tion (RPA), which are known to select only one type of 
instability. The bosonization approach and exact solu- 
tions based on the Bethe ansatz are extremely powerful 
in the purely ID case, but arc of limited use in higher 
dimension where FL quasi-particles become relevant ex- 
citations. The renormalization group (RG) technique for 
interacting electrons, however, can deal in a controlled 
way with both TLL and FL behaviors in weak coupling. 

In this respect, on the basis of the Kadanoff- Wilson 
RG scheme, Duprat and Bourbonnais (DB) applied a 
two-variable momentum approximation for the coupling 
constants, from which the competition between density 
wave and superconducting states in Q-ID systems can 
be analyzed.^ By calculating the instability of the nor- 
mal state towards an ordered state below the crossover 
scale Tx, they obtained the emergence of anisotropic su- 
perconductivity from spin-density-wave correlations due 
to the presence of nesting deviations. These effects were 
not treated in previous approaches, which are based on 
a sharp two cutoffs scaling procedure.^ 

At the end of the seventies, Lee et al. discussed the 
role of interchain (Coulomb) coupling using the RG tech- 
nique but where t±_ was excluded.^ In their work, not 
only density-wave and superconductive response func- 
tions were discussed, but also the uniform response, such 
as the compressibility and magnetic susceptibility. When 
the interchain hopping is present, however, further im- 
provement of the RG method is needed to describe fluc- 
tuations that contribute to the uniform responses like 
magnetic susceptibility, compressibility and the Drude 
weight. 

It is the purpose of the present study to further develop 
the Kadanoff- Wilson approach in the case of the Q-ID 
interacting electron systems to retain the full transverse 
momentum- dependence for the flow of coupling constants 
— a scheme that shall be denoted a A'^-chain RG in the 
following. We shall concentrate on the RG calculation of 
the uniform magnetic susceptibility and investigate the 
dimensional crossover behavior for this physical quantity 
that is directly relevant to experiment. The temperature 
dependence of the magnetic susceptibility x(T) is calcu- 
lated to show how x(^) varies with the interchain hop- 
ping t^, long-range interactions and nesting deviations 
which are well known to be relevant in organic conduc- 
tors. 

The paper is organized as follows. In §2, the model 
Hamiltonian is introduced with the formulation of the 
A^-chain RG method. The RG flow trajectories for the 
couplings constants across the crossover region arc given 
in §3. In §4, the x^) behavior in the Q-ID systems is cal- 
culated within two different approaches. We flrst consider 
the RPA method, which gives x{T) to lowest order in t±. 
The results are compared to the full A^-chain RG treat- 
ment, where t±_ is treated as non perturbatively. In §5, 
the results are discussed in connection with experimental 
data of (TMTCF)2Ar compounds and we conclude in §6. 
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2. N-chain Renormalization Group 

2.1 Model 

We consider a linear array of A^ weakly coupled chains 
of length L at quarter filling, whose Hamiltonian H = 
Ho -\- Hi consists of a non interacting and interacting 
parts 

Ho = Yl [-*4,j,aCi+l,j,<T - tj_clj^^Cij+i,^ -\- h.c. , (1) 

-ff I = ^ [Uuij^-^nij^l -\- ViHijUi+ij -h V2nijni+2j] ■ 

(2) 

Here Cij.a {c] j g-) is an annihilation (a creation) operator 
of an electron on the i-th site of the j-th chain with spin 
1 1)7 and t {t±) is the intrachain (interchain) hop- 
ping. The intrachain interaction part is of the extended 
Hubbard form with U as the onsite Coulomb repulsion, 
and Vi (V2) as the nearest- (next-nearest-) neighbor re- 
pulsion. The density operators are Uij^^ = and 
iT'ij = ^i.j.T+^^i-i.i- "^^^ kinetic term Hq can be expressed 
in the weak-coupling region as follows: 

Ho=Y^ ^p(k)Cp,k,<TCp,k,a, (3) 

p,k,(T 

U^)c^VF{pk-kF)+U{k±), (4) 

(fc_L ) = — 2f_L cos k± — 2t±2 COS 2k±, (5) 

where vp '^s the Fermi velocity and the lattice constants 
are set to be unity. The operator Cp,k,o- {Cp^^) anni- 
hilates (creates) an electron with spin a close to the 
Fermi surface of the right k ~ +kF {p = +1) and left 
A: ~ — /ci? (p = — 1) branches. The second-harmonic term 
in eq. (5), which denotes the nesting deviation, is given 
by t_L2 = t\coskF/^ts,vc? kp within the tight-binding 
approximation, however, we consider t±_ and t±_2 as inde- 
pendent parameters in this paper in order to distinguish 
the problem of the dimensionality and the nesting devi- 
ations. As it will be shown later, t±_2 does not affect the 
magnetic susceptibility. 

The partition function Z{= Tre~'^^) is written in the 
functional integral form 

Z = jj ViP*Vi,e^^-*' (6) 

where the V's are the Grassmann fields. The action S is 
given by 

s[rM = So[rM+Si[rM, (7) 
^o[V'*,V'] = E [®p(^)]"V;.(fc)V'p,.(fc), (8) 

p,fc,cr 

Si[i)*,i>] = -j^^[gpPp{q)P-p{-^) 

p-A 

+ 9aSp{q) ■ S-p{-q) 

+ 9ippp{q)pp{-q) + giaS;{q)s;{-q)\ , (9) 

where k (q) indicates the vector in the Fourier-Matsubara 
space k = (k, a;„) {q = (q, Wm)) with uin = (2n -|- l)7rT 
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{uj,fn = 2m7rr). Here, the charge and spin-density opera- 
tors of the branch p are given by 

Ppi^) ^lY.^;,A'' + q)M~k), (10) 
Spiq) ^ ^ E ^;.«(^ + q)'^"^M~k), (11) 

k,a,l3 

where ipp^a{k) is the Grassmann field corresponding to 
Cp,k,CT- The bare Green's function is given by &p{k) — 
[iujn — 'Cp(k)] ^- The couphng constants at quarter filhng 
[kp — 7r/4, vp ~ V^t) are expressed in terms of U, V\ and 
Vi as follows 



= + 4^1+61^2, 

g4p = t/ + 4^1 +41^2, 
g^a = -U. 



(12) 

(13) 
(14) 
(15) 



2.2 Renormalization Group 

In the present A^-chain RG, we achieved improvements 
for the Kadanoff- Wilson RG method by DB7 The mam 
improvements are made on the following points. 

(1) The non-logarithmic channels, which are the 
particle-hole (Landau) and the particle-particle 
(Cooper-|-) channel on the same branch, are taken 
into account. The non- logarithmic channels have been 
neglected in the previous RG, but become important 
at finite temperatures, since the logarithmic terms are 
suppressed by the temperature and become comparable 
to the non-logarithmic terms. These non-logarithmic 
terms contribute to the renormalization of g^p^iai and 
are in turn involved in the calculation of magnetic 
susceptibility. 

(2) The couplings depends on full (three) momenta 
(fci^, fc2_L, A:3j_); the forth momentum fc4j_ is determined 
uniquely by the momentum conservation. Here, ki and 
k2 are the incoming momenta and ka and k4 = ki + 
k2 — ks are the outgoing ones. When we consider the 
only fixed ki = k2 for the Cooper channel and ne- 
glect the non-logarithmic channels, two variables ki and 
ka are enough to describe the momentum dependence 
of the couplings. This scheme has been employed by 
DB. However, when we consider all scattering channels, 
namely the Cooper, Peierls, Landau and Cooper-|- chan- 
nels, the couplings can be successfully described by three 
(ki, k2, ks) instead of two variables . In order to investi- 
gate the temperature dependence of the physical quanti- 
ties, we calculate the RG at the finite temperature where 
the Landau channel plays an important role. So we keep 
three variables in the present Ai"-chain RG technique de- 
spite the complicated formulation. 

Based on the Kadanoff- Wilson RG procedure, we cal- 
culate the effective action S[Tp* ^ij^l^ at the renormalized 
bandwidth Eo{£). The partial integration is carried out 
for momentum located in the energy shell Eo{£)d£ on 
both sides of the Fermi level (Fig. 1), where Eq{£) = 
Eqc^^ is the renormalized bandwidth at step £ of inte- 
gration with the initial bandwidth Eq = 2vpkp ~ 2t. 
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Fig. 1. Illustration of the partial integration. The initial states 
(1 & 1') are scattered to the excited states (2 & 2'), and then 
to the final state (3 & 3'). The contribution of the excited state 
located in the energy shell Eq{1)M/2 are renormalized in the 
"new" renormalized action Sltj)* ,il)\ij^dl- 



By assuming the partition function invariant under the 
RG procedure, the renormalized action is obtained from 

Z = e'^(^) / / Vil}*Vil:e^^'l'' ^'l'^< [ [ V4)*Vi)e^^'^' 



(16) 



where A{£) corresponds to the free energy density at the 
step £. The partial integration with respect to if) and ^* 
can be performed by making use of the linked cluster 
theorem. Then the renormalized action is expressed as 



n=l ' \ ^i=l ' I 



(17) 



(18) 



where S\^i has i = 1 4 t/j's in the outer shell, and 

_ //PVi*P^(...)e^°[^''-^l 
^' ■ " JJV^l}*V^eSo[i'',^] ■ 

The renormalizations of gp.a,ipA(7 at the one-loop level 
{n = 2) are obtained from the outer-shell averages of the 
quadratic term of the action, i.e., ^((5*1.2)^). Depending 
on the way of choosing two among four, 5i^2 is decom- 
posed into four components as 



Sl,2 — S12 + S12 + S12 + 



(19) 



The corresponding diagrams are shown in Fig. 2. The 
definitions of are given below. Each S12 is orthog- 
onal, so that the contractions in each channel /i can be 
considered independently. Finally, the renormalized ac- 
tion is obtained as 



s[r,4>]i+di ~ s[r,i^]i = \ E((^!!2)'>, 



(20) 



where /i = C, P, L, C+ are related to the Cooper, Peierls, 
Landau, and Cooper+ channels, respectively. 

The renormalized action at the one-loop level is ob- 
tained from the shell averages of the quadratic term of 



(a) S,f (b) S, / 

k| ^ ^ ki ki > ki kr 
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+ g4±{ki±.,q'± - ki±, k3±)Alj_^p{q', k3)A4±,p{q', ki) 

+ c.c] , (24) 
-k2'' ^-kj qc-ki^' 'fqc-k3 kj-Q/' k,-Q where 

different branch hooper Peierls _ 

(,c=k,-k,) (Q=k,.k3) A,,^j,{qc,k) = Y,{i<^v)^-pA-k + qc)Mk), (25) 



(C) 

ki-j >ki+q 



Q,/3 



Landau 
(q=k3-k,) 



Ats,p(<7c, fc) = ^{io-o-y)tp-p^i3(-k + qc)ijp^a(k), (26) 

Q,/3 

a,p(g, fc) = 51 ^-p.-(^ + (27) 



Si (d) S,/ (e) S, ,c+ 




o,,p(g, fc) = ^ Cp.Jfc + 0)<T"'5Vp./3(^), (28) 

a,/3 



k,/^ ^ k4 k2-f >■ k2-q q'-k]/^ >-q'-k3 

Landau Cooper+ ^ 

(q=k3-k,) (,'=k,.k2) Opjq,k) = -Y,i'*pA'' + q)i'pAk), (29) 



same branch 



Fig. 2. Diagrams of 5i,2 for (a) Cooper, (b) Peierls, (c) Lan- ( ~ l\ = ^ \^ I * (h M ~\ °'0 I (i\ Cin^ 

dau representations on the different branch, and (d) Landau, (e) '^'y,p{q, K) — ^ 2_^y)p^^{K + q)(T tpp^i3[K), [6U) 

Cooper+ representations on the same branch. The solid (dashed) Q!,/3 

line indicates right (loft)-going electrons. The slashed lines refer ^ 

to a Grassmann field in the outer shell, A4||_p(g', k) ^ ^ ^ ''''^'^ ''Pp,a{q' — k)'4'p,aik), (31) 

(7 

the action 6*1,2, which is given in eq. (19). The contrac- ^4±A^'^ ^) ^ ^ X! '^'"''^p-'^i^' ~ ^)V'p,<t(<?')- (32) 
tion is performed by the outer shell average of all pairs 

of two V''s. Each one is given by Here qc = ki - k2, Q = k2 + ka, q = ka - ki and 



- - _ 

Si 2 — ^ 5Z \9ss{ki±: ki± — qc, k3±) In terms of eqs. (12)-(15), the above couphng constants 

icis.qc P are written as 



q' = ki + k: 



2- 



',P 

+ gts{ki±, ki± - qc,k3±)A;^ p{qc,k3)Ats,p{^c, ki) ( g^{ki 2 3±) 

+ C.C.], (21) \9s{kia,zi_) 

/.94||(^l,2,3±) \ ^ 1/1 1 W54p(fcl,2,3i) 



1/-1 -3\ /gp(fci,2,3i)\ 

il-1 1 j U^(fci.2,3^)r^ 



5^2 = -— ^ 5][gc(fcl±,Q±-fc3±,fc3±) V54±(fcl,2,3±); 2VI -lA34.(A:i,2,3±)/ 



^.^.J Cooper channel 

X O* p(Q, fc3)Oc,p(<3, fci) The contraction of the Cooper channel at the one-loop 

level is carried out as 



+ .9s(fcl±,Q± - fc3±,A:3±)Cs*p(<3,fc3)Cs,p((3,fcl) 

+ C.C.], (22) 



2 



7^ E E l^E^c(gc,gc±,fcl;^) 



'5'J;2 = -;^ E Etf''('=l-L'^2±,g± + fci±) ^^,.===t=rr . .^^'^'"F u, 



y. gu{ki±,ki^ - qc± , k'^ ; tjg^, {k'j_,k'^- qc± , A:3± ; ^) 
xa;„(,^^,)C^P,-p(9,fc2)^ ^ ^ xA:(,~c,fc3.)A.(9c,fct.), (36) 

+ ga{ki±,k2±,q± + kij_)0* p{q, ki)Oa-p{q, where 

i?o(^) 



+ 54p(fci±,fc2±,9± + fci±)C;,p(g,A:i)Op,p(g,A:2) /(.(gc, 9c_l, /ti; ^) = V -— — — — - 

4 ^ £;o(£) + AAc(gc,gc±,fcl) 

+ g^Aki^, fc2±, 9± + fci±)o::p(g, kM^piq, k2) 

T ,,,, /, I Eo{£) ^ Eoi£)/2 + XAciqc,qc±,k'A 
+ c.c] , (23) X I tanh + tanh — ^ 

y ^ ^ 



(37) 



LN 

~ki,a,q' P 



^c{qc,qc±,k'A = vpqc +C±(fc±) -C±(9c± - fc±), 

(38) 

XA411 (g',fc3)A4|| (g, fci) N c n ^ ^qan 

"'^^ " Vp'O'C = ^±(fci_L - gc_L) - ^±(fci±)- (39) 
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Fig. 3. Plots of /c(qc = 0), /p(Q = Qo) (solid line), and /L(q = 
0), Ic+{<i' = <?()) (dashed line) as a function of Eo(^)/4T, where 
Qo = {2kp,iT) and qg = (2fcj?,7r). When the fi channel exhibits 
complete logarithmic divergence, 1^ reaches unity. 



The function Iq indicates the degrees of logarithmic di- 
vergence; the case Ic = 1 imphes a Cooper instabihty 
showing a complete logarithmic divergence, while it is 
partially suppressed for /c < 1 (Fig. 3). This Cooper 
channel is shown diagrammatically in Fig. 4(a). We note 
that qc is set to be zero in DB's approach assuming that 
the Cooper instability is always logarithmic. 

2.2.2 Peierls channel 

At the one-loop level [Fig. 4(b)], the contraction in the 
Peierls channel is given by 

1 

2^ 



T 

LN 



E E 



de 



Nttvf 



^/p(0,Q±,A:l;£) 



xOt{Q,h±)0,{Q,k±), (40) 
where 



E 



1 



X tanh — — — h tanh 
4T 



^^^^Eoie) + XAp{Q,Q^,k'^) 

Eo{e)/2 + XAp{Q,Q^,k'J' 
2T 

'(41) 



ApiQ,Q^,k'^^) = vfQ + C^ik'J + C^ik'^ - Qi_), (42) 
VFQ^-^±{k3±)~UiQ±-k3^). (43) 
The flow of /p is similar to Ic (Fig. 3). 

2.2.3 Landau channel 

The one-loop contraction of the Landau channel gives 
a non- logarithmic contribution to the RG equations. This 
channel was discarded in the previous RG. On the tem- 
perature scale, however, the non-logarithmic terms be- 
come effective when the Cooper and Peierls logarith- 
mic singularities are suppressed. Moreover, the spin- and 



charge-coupling terms (e.g., 5pg4cr, <74pff4cr) give contribu- 
tions of the same order of the spin- and charge-decoupling 
terms (e.g., g^gaQia) 

The Landau contraction [Fig. 4(c)] yields (see Ap- 
pendix A for the details) 

1 

2* 



T 



E E 



di 



Y,lL{q,q±,k'^;£) 



LN ^ ^ Nttvf ^ 

{fel,2:9}* 

X .g;y (fcl± , fcl - 9_L , fcl± + 9_L ; ^)5i/' (fel , k2± ,k'j_- q±;£) 

xO:{q,h±)0,,{q,k2±), (44) 
where 

1 



A tanh 



.A=±l 



4 ALiq.q±,k'^) 

Eoi£)/2 + XAi^iq,q^,k^) 
2T 



(45) 



Ai^iq, q^,k'^)^ VFq + U{k/^) - ^±{k'^ -q±), (46) 

VFq = U{ki±)-S.±{ki± + q±). (47) 

The function /l does not reach the unity, i.e., the Landau 
channel is always non logarithmic (Fig. 3). 

Similarly, the contraction of the Cooper-|- channel at 
the one-loop level [Fig. 4(d)] is 



di 



— y y 

T.N ^ ^ 

i/,z.'=4||,4±{fe^3,,~'}. 



LN ^ NnvF ^ 



Y,Ic+iq',q'±,k'^;£) 



X gi,{kij_,q'j_- kij_,k'j_; £)g^' {k'^, q'^ - k'^, k^^ ; i) 

X^l[q'M±)^u'{q'M±), (48) 
where 

Eo{£) 



Ic+iq',q'^,k'^;£) 



1 



, (49) 



4 Ac+{q',q'^,k'^) 

^^ ^^ Eo{£)/2 + XAc+{q',q'^,k'J 
2T 

LA=±l 

Ac+iq\q'^,k'^)=VFq' + ^Akl)+^±{k'^-q'^), (50) 
VFq' ^ -Uiki±) - U{q'± - ki±). (51) 

2.2.4 Flow equations 

Taking into account all the above one-loop contrac- 
tions (Fig. 4), we can finally write the flow equations in 
the form 

■^G^{ki±,k2±,k3j_) = ^y^Ft,{ki±,k2±,k3j_,k'^). 



(52) 
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(a) Cooper 
k, 



(b) Peierls 
k. 



k| -Q 



k' Q 



k, Q 



(c) Landau 
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- k'+ q 


. -k2 
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(d) Cooper- 

-k, +q' 



- k" + q' -kj + q' 



Fig. 4. Diagrams considered in the flow equations at the one- 
loop level. Here, qc = ki — k2, Q = k2 + ka, q = ka — ki and 
q' = ki + k2 



Here, v denotes the type of the couphng, v — p,a, 4p, 4ct. 
The function F^, contains the quadratic terms of the cou- 
phng Gu{— gvl'KVp) and the contributions of the loop 
The F,y's are given by 

X [Gp (fcii, /cij_ - qc^,k'^)Gp{k'^,k'^ - gc±, ^3^) 

+ ^/p(Q,Q±,fc±) 

X [Gp{ki^,Q±_ - k'^,k'^)Gpik'j_,Qj_ - ks^^k^^) 
+ 3G^{ki^, Qi_ - fcl, fcl)G,(fcl, gj. - fc3±, fc3±)] 

X [-Gp(fci±, fc^ - qj_,qi_ + ki±_)G4pik'^,k2±,k/j_ - q±) 

- G4p(fci±, fc^ - q±,q± + ki±_)Gp{k'^, k2±,k'j_ +q±;e) 
+ Gp{ki±,k'^ - q±,q± + kn_)Gi„{k'^, k2±, k'^ - q±_) 

+ G4cr(fci±, fc± - q±,q± + ki±)Gp{k'^,k2±,k'j_ + q±)] , 

(53) 

Fcr{ki±, k2±,k3^, k'^) = i/c(9c, 9ci, ^1) 

X [2Gcr(fci^,/si_L - gc±, fcl)Gcr(fcl, fcl - qc±,ks,A_) 

- Gp{kij_,kij_ - qc±.,k'j_)Gaik'j_,k'j_ - gc±-fc3i) 

- Gcr{kij_,kij_ - qc±, k'j_)Gpik'j_, k'^ - qc±, k3±)] 

+ \lp{Q,Q±-k'^) 

X [2G,(fcii, - fcl, fcl)G,(fcl, Qi_ - fc3±, A:3±) 
+ Gp(A:i_L,Q± -fcl,fc^)G^(fcl,Oi -/c3_l,A:3_l) 
-I- Gct(A:i_l, (5_L - k'j_, k'j_)Gpik'j_, Q± - k^^.k^^)] 



+ ■^iL{q,q±,k'j^) 

X [Gaiki±,k'j_ - q±,q± + kij_)G4p{k'^, k2±, k'^ ~ q±_) 
+ G4p(fci_L, k'j_ - q±,q± + ki±)Ga{k'j_,k2±, k'^ + q±;£) 

- Ga{ki±,k'j^ - q±,q± + kij_)G4a{k'j_, k2±, k'j_ - q±) 

- G4a{ki±, k'^ - q±,q± + ki±)Gaik'j_, k2±, k'j_ + q±)] , 

(54) 

F4p{ki±,k2±, k3±,k'^) = 2^L{q, q±, fcl) 

X [-2Gp{kij_,k'^ - qj_,qj_ + ki^)Gp{k'^,k2j_, k'j_ - q^) 
+ G4p(fci^,/j^ -qj_,qj_ + ki4_)Gip{k'^,k2±.k'^ + q±;i) 
+ G4cr(fci±,fcl - q±,q± + kij_)G4:aik'j_,k2±, k'j_ - q±) 
+ G4,p{ki±,k'^ -qj_,qj_ + kij_)G4:aik'j_,k2±,k'j_ + q±) 

+ Giaiki^^k'^ - qj_, q± + /sij_)G4p(fc^, fej., k'j_ + q±)] 
1 

~ T 

X [G4p(fci 

+ G4.(fcl±, -z^ - k[^,k'^)G4a(k'^,q'l_ - k'^,k3±)] , 



-Ic+{q ,q'j_,k'j^) 

- fci±>fcl)G'4p(fcl,g^ - k'^,k3±) 



(55) 



F4a{ki±,k2±,k3±,k'_^) = ^/l((7, (7_L,fc^) 

X [~3Ga{ki±,k'^ - q±,q± + ki±)Ga{k'j_, k2±,k'j^ - qj_) 
+ G4p{ki±,k'^ - q±,q± + ki±)G4„{k'^, k2±, k'^ + g^) 
+ G4cr(fci±,fcl -q±.q± + ki±)G4p{k'^,k2±,k'j^ + q±)] 

-i/c+(g',gl,fc^) 

X [G4p(fci^, 9^ - k[j_, k'j_)G4<j{k'^,q'i_ - k'^, fcax) 

+ G4^(fci±, - k[^,k'^)G4p{k'^,q'^ - fc^L, ^3±)] , 

(56) 

If we neglect the fcx-dependence, these flow equations 
agree with those obtained in the ID limit [eqs. (11)-(15) 
in Ref. 10]. Calculating simultaneously these differential 
equations, we can obtain the renormalized couplings at 
each temperature. For the calculation of the A^-chain sys- 
tem problem (2A^ points on the Fermi surface), we solve 
4A^^ differential equations. 

3. Flow of couplings & crossover behavior 

In this section, we show the results of the numerical 
solutions of the flow equations (52) ^ (56) for N = 31- 
chain systems with U = 3t, Vi = V2 = 0, t± = O.Olt 
at T = 10~H, and analyze the flow of the couplings. We 
set the minimum of the renormalized bandwidth, Eq{£c), 
of the order of 10^^i?o, which is low enough for T = 
10~^ ~ 10~^t. For N >9, the results become essentially 
independent of the number of chains. The present N- 
chain RG exhibits a clear difference in the property of 
spin gap between the odd-number chains (gapless) and 
the even- number chains (gapped) . The aim of the present 
paper is to investigate the many-chain systems, where 
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Fig. 5. Illustrations of the coupling C, P, C+P, L. 



(a) 




le-05 



le-04 0.001 0.01 



(b) 




le-05 le-04 0.001 0.01 



Fig. 6. (Color online) Energy flow of the couplings Gp (a) and 
-Ga (b) for U = 3.0t, Vi = V2 = 0, = O.Olt at T = 10""'*. 



the spin-gap does not open, so that we focus on the odd- 
number-chain systems in the present paper. 

The renormahzed couphngs depends on three vari- 
ables {ki±, k2±, k3±). In order to see the essential points 
of the renormalized couplings, we show the energy flow 
of the couplings for some typical pairs: 

(C) Gi/(0, 0, 0): the scattering of the Cooper pair ki± = 
k2± = 0. 

(P) Gu{Tr, 0, tt): the scattering of the Peierls pair ki± ~ 
TT and —k2i_ — 0. 



(C+P) G,.(7r/2,7r/2,7r/2): the scattering of the pair 
fci_L ~ 7r/2 and — fc2± — ~7r/2, where the Cooper 
and Peierls instabilities coexist. 

(L) Giy(7r/2, 0, 7r/2): the scattering of the pair fcix — 7r/2 
and —k2i_ — 0, where both the Cooper and Peierls 
instability are weakened. 

Each pair states are illustrated in Fig. 5. Note that, 
in the actual calculations, tt and 7r/2 in the arguments 
of the above coupling constant are replaced as tt ^ 
(30/31)7r, 7r/2 — > (16/31)7r due to the finite number of 
chains. 

The flows of the corresponding couplings arc shown 
in Fig. 6. For i?o(€) > At±, each couplings are indistin- 
guishable and fully agree with the couplings of the ID 
case (the dash-dotted line in Fig. 6). On the other hand, 
for EQ{t) < At_i, a difference becomes noticeable, indi- 
cating the behavior of the interaction is no more ID, but 
Q-ID. Thus, we can flnd the characteristic energy of the 
crossover as = 4t_L. This is a clear evidence that the 
structure of the interaction in Q-ID correlated electrons 
do exhibit the crossover behavior. 

Here a few remarks are in order about the crossover 
behavior and temperature appearing in the response func- 
tions. In the presence of t_L, there is a characteristic 
temperature for response functions, which are calculated 
through the RG equations. With decreasing temperature, 
the effect of t± on the RG equations becomes relevant, 
suggesting a change of the state from the ID TLL regime 
into the 2D regime. Generally, t±^ in the presence of in- 
teraction is renormalized to a reduced value, and the 
crossover temperature is also renormalized as T* ~ ty_'^ 
as obtained from the two-loop RG with the self-energy 
corrections.^' Here we distinguish T* {E*) from the 
bare one. There is no renormalization for Tx('^ t±) at 
the one-loop level, i.e., no self-energy corrections are in- 
cluded. Nevertheless, a crossover behavior does exist in 
the two-particle quantities at this level of approximation 
as shown in Fig. 6. The effect of t± on the two-particle 
correlations such as those contributing to the magnetic 
susceptibility can be thus obtained even at the one-loop 
level and the results are expected to carry over to the 
case where Tx is replaced by T*. 

In the Q-ID regime, where Eq{£) < E^, the present N- 
chain RG shows that the coupling C (P) is reduced (en- 
hanced) compared with the ID situation. In ID systems, 
(fc, —k) pairing contributes to both Cooper and Peierls 
instabilities. In the Q-ID case, however, the Cooper pair- 
ing such as in C of Fig. 5 is characterized by bad nesting 
conditions and contributes weakly to the Peierls instabil- 
ity. The C pairing thus flows to the attractive sector and 
is attractively relevant. Similarly, for some Peierls pairing 
such as in P, the contribution to the Cooper instability 
is weak so that the coupling P flows to the repulsive sec- 
tor and is repulsively relevant. On the other hand, for 
Cooper pairs like the C+P, both the Cooper and Peierls 
instabilities arc involved even in the Q-ID regime. In this 
case, the coupling C+P becomes marginal below E^. Fi- 
nally for the pair L, both Cooper and Peierls instabilities 
are weakened so that coupling L is hardly renormalized. 
Note that even for the coupling L or C+P, the couplings 
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(a) Gp {kfj_, k2j_, k,J 




-n 71 



Fig. 7. Renormalized couplings Gp{ki±,k2±,ki±) (the upper 
panel) and —Ga-{ki±, k2±,ki±) (the lower panel) for U/t = 
3.0, Vi/t = V2/t = 0.0, tj_ = O.Olt at T = 10"**. 

are slightly renormalized in the thermal shell Eo{£) ~ T 
due to the non-logarithmic terms, i.e., the Landau and 
Cooper+ channels. (See Eo{£)/Eo - lO^* in Fig. 6.) 

In Fig. 7, we plot the renormalized couplings 
G,y(fcij^, fc2i, /cij^) along the plane ki± = k^^ 
in the parameter-space A:2_l, ^3^)- The function 

Gv{kn_, k2±, ki±) is the coupling in the limit of q ^ 0, 
and thus corresponds to the Landau's interaction func- 
tion /(ki,k2), which is the second functional derivative 
of the energy per unit volume.^ As we can see from Fig. 
7, there are two remarkable structures. One is the valley 
along ki± — k2± (the valley C), indicating the Cooper 
instability of which momenta are ki and — k2. Another 
is the ridge along fcix = —k2i_ ± tt (the ridge P), indi- 
cating the Peierls instability reflecting the nesting vector 
Qo = (±2fcF,±7r). If the system is pure ID {t± = 0), 
the value of the renormalized couplings are (Gj^)* ~ 0.7 
and {-G^J^y ~ 0.1. The valley C is smaller and the 
ridge P is larger than {G^P)*, indicating that the valley 
C is more attractive and the ridge P being more repul- 
sive. Moreover, the height of the ridge is not uniform. 
The cosine-like curve of the valley C suggests the d- or 
/-wave instability, and the maximum (minimum) of the 
ridge P suggests the hot (cold) spots of the density-wave 
instability.^ Namely, we can obtain the information of 
the symmetry of the Cooper pair and the structure of 
the density-wave fluctuation only from the renormalized 



couplings. For the quantitative comparison of each in- 
stabilities, we need to calculate the response functions, 
which are presented in Appendix C 

Interestingly, the structure of the present couplings 
in Q-ID is quite similar to the /(ki,k2) for the two- 
dimensional (2D) Hubbard model near half-filling.^^ This 
suggests that the electronic properties of the Q-ID sys- 
tems have a lot of similarities with those of 2D systems 
near half-filling. For example, the Landau parameter Fq 
can also become positive in the Q-ID case, which sug- 
gests the possible existence of the spin-dependent zero- 
sound (zero-spin-sound) . This possibility is supported by 
the decrease of the magnetic susceptibility at low tem- 
peratures as shown in the next section. Furthermore, it 
might be possible that becomes negative, which corre- 
sponds to the reduction of the Drude weight in the Q-ID 
case. 

The possible fixed points of the couplings in the param- 
eter space {ki±, k2±, fcsj^) are summarized in Fig. 8. The 
couplings on the C plane kn_ — k2±, are attractively rele- 
vant due to the Cooper instability. On the P plane where 
k2i_ + fcsx = iTT, the couplings are repulsively relevant 
due to the Peierls instability. At the intersection of the C 
and P planes leading to the C-l-P line, the couplings are 
marginal below due to the coexistence of the Cooper 
and Peierls instabilities. Remaining couplings are hardly 
renormalized in the Q-ID regime. These conditions for 
the coupling constants have to be taken into account in 
order to construct an effective low-energy Hamiltonian 
in the Q-ID regime. 



plane C 
(attractively relevant) 




(repulsively relevant) 

Fig. 8. (Color online) Illustration of the total structure of the 
fixed points for the couplings Gp and Ga in the parameter space 
(fcixi t^2±) k^^). The couplings on the plane C (fcix = fc2±) a^r^ 
attractively relevant due to the Cooper instability. The couplings 
on the plane P {k2± + = ztw) are repulsively relevant due 
to the Peierls instability. Along the line C-l-P, the couplings are 
marginal below Tx due to the interference between Cooper and 
Peierls instability. 
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4. Magnetic susceptibility 

In this section, we investigate the properties of the 
magnetic susceptibihty. We first briefly review the prop- 
erties of x(r) in the ID hmit^" and then tackle the main 
problem of x{T) for Q-ID systems by first considering 
the lowest order perturbation effect of t± with respect 
to the ID RPA solution for {t±/T)'^ < 1, and second by 
the full one-loop RG treatment where t± is treated non 
perturbatively. 

4-1 Magnetic susceptibility in ID 



(a) 



s 

o 





V\ increases 




j/=4r 
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-^V\ V2=0-5-2r 




--. ^^^^ ■■^Svvv^x 




~~~~~~ "^,;^^VVVVv 


increases "^^-^^^^^N^^^$^^ 







0.2 0.4 0.6 

Tit 




Fig. 9. (Color online) Magnetic susceptibility of ID, x^^(T), as a 
function of temperature, (a) The black-solid line denotes (T) 
for U = At, the gray-solid lines denote that for U = 4t, Vi = 
0.5 ~ 2t, the gray-dashed lines denote that for U = 4t, V2 = 
0.5 ~ 2t. (b) Comparison of x^°(T) by RG to the other method. 
^ID ^rp^ RPA is shown by dashed-lines. The maximum values 
of X^^C^) by QMC locate around the dash-dotted-lines.i^ The 
exact solution at T = are indicated by the square and circle. 



We plot the ID magnetic susceptibility in Fig. 9 
(a) for iU,Vi,V2) = (4t,0,0) (thick solid-line), (4^,0.5 - 
2t,0) (thin solid-lines), (4t, 0,0.5 - 2t) (thin dashed- 
lines). (The results coincide with those of Ref. 10, or 
those of the following section in the limit t± =0.) For 
the Hubbard model {U >0,Vi = V2= 0), x^^iT) is pro- 
portional to 1/T at high temperatures, has a maximum 
around T ~ 0.2t, and then decreases towards the zero 
temperature value x^°(0) = i2/7rvF) [I - U/{27rvF)]~^ 



with an infinite slope. When Vi increases, X^^{T) is en- 
hanced but the temperature profile remains essentially 
the same; the influence of Vi being similar to the one 
of U. When V2 increases, however, x^^{T) decreases and 
becomes flat at low temperature, that is dx^^ /dTlr^a = 
for [/ = 2V2 (i.e., ga = 0)- 

In Fig. 9 (b), we compare the temperature profiles of 
X^^{T) for different approaches. The results of the RG 
are shown by the solid lines and the ID (renormalized) 
RPA form for x.^^(T), 

2 Xo{T) 



(57) 



TTVp 1 + [G, (T) - G4± (T)] xo {T)/2 ' 

are given by the dashed lines. Here xo (T) = tanh(i?o / 4r) 
is the bare magnetic susceptibility per branch p = ± 
normalized by ttvf- Ga{T) denotes the renormalized 
coupling at T, which is given by Ga{T) = GafW -I- 
G<,ln(2r/i;o)], and G4_l = (G4p - G4^)/2. We can see 
that eq. (57) includes the fluctuation beyond the ordinary 
RPA, where Ga{T) = G^ is independent of tempera- 
ture^^ . The maximum values calculated by the quantum 
Monte Carlo (QMC) technique^'^ (for which low temper- 
atures cannot be reached due to sign problem) are in- 
dicated by the dash-dotted lines. The exact solutions at 
T = are shown by the square and the circle for U = 2t 
and At, respectively.^'^' 

For U = 2t, x^°(T) of the RG and renormahzed RPA 
results, eq. (57), quantitatively agrees with that of the ex- 
act solution and QMC. However, for U = At, Xrpa(-^) 
sizably larger than the QMC, suggesting that RPA result 
overestimates the effect of spin correlations. This overes- 
timation is not found in the full one-loop RG result, and 
X^^(T) remains close to that of QMC. This improvement 
is due to non- logarithmic terms, which becomes effective 
at high temperatures, and also comes from the mode- 
mode coupling terms. 

Summarizing, the RG prediction for x^'^(T) agrees 
well with the exact solution at T = even for large U. It 
also agrees reasonably well with the QMC results at high 
temperatures, though it only slightly departs from QMC 
at large U. Therefore the present RG technique yields 
the most reliable results for the uniform susceptibility in 
the case of the ID extended Hubbard model for the wide 
temperature region < T < i. 

4-2 RPA approach with the perturbation oft± 

First we see the variation of the couplings Gi, on 
the basis of the perturbation of t±. We adopt the ID 
Kadanoff- Wilson RG of Ref. 10, as a non-perturbative 
part. The functions /c,p.l,c+ can be expanded in t±_ (or 
A^) as 



tanh 



AT 



T 1,-2 

Il,c+ - -^jT- cosh 



-fiEoie)/AT) 



A, 



C,P 



2T 



(58) 



AT 



9iEo{e)/AT) 



A 



L,C+ 

2T 
(59) 
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where 







1 + 









(a) magnetic susceptiblity 
,k 



x2 — cosh(2a::) 
3 cosh^ X 



(60) 
(61) 



d/ 



k+qo 



k 

+ 

k+qo 



k'+qo 



k+qo 



(Here, we omit the (g, <7_l, A:j_) dependence of and A^.) 
The flow equation (52) for q — > 0, which is relevant for 
the magnetic susceptibility, can be rewritten as follows. 




k'-Qo 



GAT) =Gf{T)-2GlHT) ("^ 



G4p(T) =GlD(r) + -{Gl^ + GL)g(r) 



G4.(T) =Gl°(r) + G4pG4.g(r) ( ^ 



T 



and for G4^ = (G4p — Gicr)/2, 

Gi^{T) ^G'R{T) + Gl^g{T) 



where 

Gy{T) = Gt,{kij_,k2j_,kij_;£,T) 



ki ,A;2 



(62) 
(63) 

(64) 

(65) 

(66) 
(67) 



(c) superconductivity 
k 




and /(T) = Jd£ f{Eo{l)/4T), g{T) = jdl g{Eo{e)/AT). 
For T ^ 0, /(T) ^ 7C(3)/47r2 ~ 0.213, g(r) ^ 0. The 
flow equations with full-momenta dependence are derived 
in Appendix B. Here we use the relation 

ki ,/c2 k' ki ,k2 k' 

= EE{^c+(<z',9l,fc^)}' 

k\ , k2 k' 



= mi 



(68) 



Therefore the enhancement of Gi, , aside from Gp , is pro- 
portional to {ti_/Tf. 

Next, we estimate the magnetic susceptibility with 
these Gj/'s. By extending the RPA form eq. (57), the 
magnetic susceptibility in Q-ID lead to be 

[xrpa(t)]-^ = [1 + [G^T) - G4^(r)} xo(r)/2] . 

(69) 

By using eqs. (63) and (66), we can write X'RPa{T) at 
{t^/Tf < 1 in the form 

[xrpa(T)]-^ ^ (xJipa)"' - ^ [AG.(T) -I- AG4^(T)] , 

(70) 



Fig. 10. Diagrams for the renormahzation of the pair vertex 
at the one-loop level. The triangles and the squares denote the 
pair vertex and the coupling for (a) the magnetic susceptibility, 
(b) the density wave, and (c) the superconductivity. The wavy 
line denotes the corresponding source-fields h^. 



Consequently, the magnetic susceptibility is enhanced by 
t±_ from that of ID system on the basis of the pertur- 
bation theory. Moreover, eq. (70) indicates that a shift 
from Xrpa exists even above Tx ~ t_L . This suggests that 
there is no clear feature signaling the crossover from ID 
to Q-ID in the temperature dependence of the magnetic 
susceptibility. This is consistent with the experimental 
results for the TMTTF and TMTSF salts.'' 

Since this is the perturbation of {t±/T), one may 
think that the perturbative form eq. (70) is not valid 
for tj_/T ~ 1, namely where the crossover behavior is 
expected to take place. Even in such case, we can es- 
timate xrpa(T) of eq. (69) from the numerical evalua- 
tion of Gv{T). As is shown later in Figs. 12 and 13, the 
Xrpa(T) so obtained qualitatively agrees with the per- 
turbative form, namely, XRPA(r) is enhanced by t± for 
the whole temperature region. 

4-.3 Renormalization group approach 

In a path integral approach, the uniform magnetic sus- 
ceptibility or the response functions can be calculated by 
adding a set of source fields /i^ to the action. The deriva- 
tion of response functions is given in Appendix C. In the 
case of the uniform magnetic susceptibility we intro- 
duce the following Zeeman coupling of a source field to 
the spin density operator: 



s':[r,i^] = y2[hp{q)o;{q,k) 



k,q 



C.C. 



(73) 



where 



The one-loop corrections yield 



AG AT) - 2Glf{T) 



T 



E hp{q)z,^p{k,k + q;£)[0*Aq,k) 



>o, 



AG4^(r)=GL5(r)(^) >0. 



(71) 
(72) 



where Za 



+ X'y{^)hp{q)hp{q). 
is the pair vertex. 



(74) 
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Fig. 11. Diagrams of the magnetic susceptibility contained in the 
Af-chain RG. The dots indicate the bare coupling constants. Di- 
agrams in the second line are considered in the simple RPA. 

For the uniform magnetic susceptibility, i.e., q ^ Qo = 
0, Zc^^p is expanded in the form^*^ 

Za.p {kj_:£+ d£) = z^^p {k±^ ; t) 

X /L(go,go±,fcl;^)dA (75) 

where 

= -G„ + (G4p - G4.)/2. (76) 

The corresponding diagram is shown in Fig. 10 (a). 
The flow equation for z^-.p is obtained as follows: 

^lnz<^,p(fc^) = ^^G^(A:±,fc^,fcx)/L('?o,9o±,fc±)- 

(77) 

The uniform spin susceptibility in units of [t:vf)^^ for 
both branches is given by 

^(^) = 1^E/ — Kp(fc±)]'^L(go,go±,fc±). (78) 

The diagrams which are included in the present 
method are shown in Fig. 11. It must be noted that 
mode-mode type of couplings (the third line of Fig. 11) 
are included in the present calculation. The contribution 
from these diagrams introduce a qualitative difference 
for x{T) with respect to RPA, as is discussed in the next 
subsection. 

4-.^ Susceptibility in the Hubbard limit 

The numerical calculation of the magnetic susceptibil- 
ity is performed by taking = 15, which is large enough 
to discuss the infinite system in present numerical result. 
First we calculate the magnetic susceptibility of Q-ID 
Hubbard model by taking tj, > 0, J7 > 0, Fi = V2 = 0, 
which is useful to understand the effect of t± for weak 
interaction. Figure 12 shows the magnetic susceptibility 
as a function of temperature for the small U{= l.Ot) 
and t± — 0.1 ^ 0.3t. The behavior of x{T) qualitatively 
agrees with the one found in ID, where there is a broad 
peak with a maximum around T 0.2t and a slight drop 
close to the zero temperature. The effect of t± is shown 
in the inset of Fig. 12, as the ratio of xiT) to x^^{T)- A 
remarkable enhancement of x(^) from x^^(^) appears 
below T - 0(t±). We find that x(T) is enhanced by 
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T/t 

Fig. 12. (Color online) Magnetic susceptibility as a function of 
the temperature, xC^): for i±/t = 0.1 ~ 0.3 (from bottom to 
top), tj_2 = with U = l.Ot, N = 15. The result by the modified 
RPA is also plotted for t± = O.lt. The inset shows the ratio of 
X(T) in Q-ID to that in ID, xiT)/x^°iT)- 
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Fig. 13. (Color online) Magnetic susceptibility as a function of 
the temperature, x(T)i for *x = O.lt (solid line), t± = 0.24 
(dashed line) and tj^ = 0.3t (dotted line) with U = 4.0t, N = 15, 
tj_2 = 0. The dash-dotted line denotes the ID case. The result 
by the RPA is also plotted for = O.lt. The maximum of Xrpa 
for t± = O.lt exceeds 4.0x0- The inset shows the ratio of x(^) 
in Q-ID to that in ID, x(r)/x^°(T). 



t± as long as U is small. Such a behavior qualitatively 
agrees with that of the RPA, though the enhancement 
of XRPA(r) is larger because the RPA overestimates the 
influence of spin fluctuation. 

However, a noticeable feature emerges for large U. Fig- 
ure 13 shows plots of x(T) and xiT)/x^^iT) for U = At, 
and t± = 0.1 '--^ 0.3t. Although the global feature is sim- 
ilar to the ID case, the qualitative difference appears at 
low temperatures. Compared with x^^{T), x{T) is sup- 
pressed at low temperatures, whereas x(^) > X^'^(^) is 
seen at high temperature. Note that this property is sim- 
ilar to that found earlier by Lee et al., who treated only 
interchain - Coulomb - coupling.^ 

Such an unexpected property definitely differs from 
the RPA result for susceptibility, which shows an en- 
hancement with respect to the ID limit over the whole 
temperature region. The difference in xiT) between the 
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present RG and the RPAQoriginates mainly from the 
effect of mode-mode couphng, which represents the re- 
sult of the quantum interference of each channel, shown 
in the third line of Fig. 11. The mode- mode coupling 
term is an higher-order interaction term, so that its ef- 
fect appears for sizable amplitude of interaction. The 
RPA overestimates the effect of spin fluctuations espe- 
cially for the large interaction. Such an overestimation 
of RPA can be corrected by the mode-mode coupling, 
as seen from the successful explanation of various exper- 
iments on itinerant electron magnetism. The present 
results for the Q-ID electron systems suggest the impor- 
tant role of the mode-mode coupling (thus the quantum 
interference of Landau and the other channels) in the A^- 
chain RG. It is worth noticing that for antiferromagnetic 
spin- fluctuations, their enhancement by the mode- mode 
coupling also leads to the reduction of x(T). This ana- 
log situation was discussed for the two-dimensional sys- 
tems with the nested Fermi surface, where xC^) is sup- 
pressed by the antiferromagnetic spin-fluctuations due to 
the mode- mode coupling. 

This behavior of Q-ID x(T) can be understood intu- 
itively as follows. At temperatures higher than t±, the 
ID fluctuation dominates and suppresses x^^{T)- How- 
ever such an effect of ID fluctuations is weakened by t±_ 
resulting in the enhancement of On the other hand, 

at temperatures lower than t± , the formation of the SDW 
(or antiferromagnetic) ordered state is expected due to 
two-particle interchain couplings. Thus such an effect of 
the high dimensionality fixing the direction of spin re- 
duces x(T) compared to the ID case. This reduction is 
specific to spin fluctuations with small momentum trans- 
fer q = 0, at variance with spin fluctuations with the 
large transfer Q = {2kp,7r) (those contributing to the 
SDW response function) which are enhanced over the 
whole temperature region by tj_ (see Fig. C-1). 

4-5 Susceptibility in the presence of long-range interac- 
tions 

We now examine the magnetic susceptibility in the 
presence of the long-range interaction Vi, which takes 
an important role in the following Q-ID organic conduc- 
tors. In the TMTTF salts for example, nearest-neighbor 
repulsion Vi is expected to be large because of the pres- 
ence of a charge ordered state. The electronic state 
of the TMTSF salts, which exhibits the coexistence of 
SDW and CDW,^^'^^ is explained by taking account the 
large next-nearest-neighbor repulsion V2-^*''^^ The spin- 
triplet superconducting state in the TMTSF salts^^ could 
be originated from the charge fluctuation, which is en- 
hanced by ]/2.^^'^^"^^ In view of these observations for 
the TMTTF and TMTSF sahs, we examine the follow- 
ing two cases. (A) f7 > > 0, V2 = for the TMTTF 
salts; {B)U>Vi>V2>0 and V2 = U/2 for the TMTSF 
salts. We choose the following parameters {U,Vi,V2) — 
(1.6t,0.7t,0) for (A) and {U,Vi,V2) = {0.8t,0M,0At) 
for (B), which give the same power law exponent of the 
SDW-response function as for {U, ^1,^2) = (4.0i, 0, 0) in 
the ID case. 

It should be noticed that the present calculation is 
performed without the Umklapp scattering, which also 
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Fig. 14. (Color online) Magnetic susceptibility as a function of 
the temperature, xC^); for *± =0.1 (solid line), = 0.2 (dashed 
line) and tx = 0.3 (dotted line) with (!7, Vi, ^^2) = (1.6t, 0.7t, 0), 
N = 15, tx2 = 0- The dash-dotted line denotes the ID case. 
The inset shows the ratio of x(^) iii Q-ID to that in ID, 
X(T)/X^°(T). 



1.28 



s 

1.16 















Tit 




0.1 0.2 0.3 0.4 X:, 








1 \ 








0.996 >i, \. 








0.992 ^ \. 
^ V 


/ "5' 


0.988 \. - 
\\ 




0.984 ^ 



0.1 0.2 0.3 0.4 



Tit 

Fig. 15. (Color online) Magnetic susceptibility as a function of 
the temperature, xC^): for *x = 0.1 (solid line), tx = 0.2 (dashed 
line) and ix = 0.3 (dotted line) with TV = 15, tx2 = and 
(C/, Vi,y2) = (0.8t, 0.5t,0.4t), where both the SDW and CDW 
instability are dominant. The dash-dotted line denotes the ID 
case. The inset shows the ratio of x(^) ™ Q-ID to that in ID, 
x(T)/xiD(T). 

exists at the quarter- filling (i.e., Sfc^-Umklapp scatter- 
ing) and may give rise to a charge ordered state. "'^^'^^ 
In ID case, the charge ordered state does not affect 
the magnetic properties due to the spin-charge separa- 
tion. Therefore the present results of x(r) without Sfc^?- 
Umklapp scattering would be valid at high temperature 
corresponding to the ID regime. The effect of Umklapp 
scattering, which is expected in QID regime, could be 
small even at low temperature, since the magnetic sus- 
ceptibility of TMTTF salts shows no clear signal around 
the charge-order transition. Thus the effect of the Sfcj?- 
Umklapp scattering is naively expected to be small in 
the present case in general. ^^'^^ However, it still remains 
as a future problem to clarify the qualitative effect of 
charge order on the magnetic susceptibility in the Q-ID 
regime. 
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Fig. 16. (Color online) Magnetic susceptibility for U = 4.0i, V\ = 
V2 = 0, A'' = 15 with = O.lt, tx2 = (perfect nesting) and 
tx2 = O.Olt (imperfect nesting). The inset shows the response 
functions as a function of the temperature. The subscripts i of 
SDW and dSS indicate the each response function for ix2 = 
{i = 0) and 4^2 = O.Olt {i = 1). 



The temperature dependence of the magnetic suscep- 
tibihty x{T) and its ratio to x^^C^) are shown in Fig. 
14 for {U,Vi,V2) = (1.6^,0.7^,0) and in Fig. 15 for 
([/, Vi, 1^2) = (0.8t, 0.5t, QAt). The SDW instabihty is the 
most dominant for the former case, while both the SDW 
and CDW instabihties become most dominant in the lat- 
ter case. 

As is clear from the plots of the ratio x{T)/x^^{T)^ 
x{T) for Vi,V2 > is also reduced from x^^i^") at low 
temperatures. In these cases (especially in the case with 
U = 2V2), the enhancement of x{T) is small, and the 
relation x(T) < x^^{T) is seen for the wide range of 
temperature, i.e. , at least for T < OAt. These features 
are understood in terms of the coupling g^i— —U + 2V2). 
Roughly speaking, the coupling go- is renormalized to 
zero as G„(T) = G^/ [1 + G^\n{2T/Eo)], reflecting the 
effect of the ID fluctuation. But for t± > 0, the renor- 
malization is stopped by t± as eq. (63), which leads to 
the enhancement of xC^) at high temperatures for the 
case of the Hubbard model with on-site interaction. On 
the other hand, for small go-, the renormalization is small 
even in the ID case, i.e., the effect of ID fluctuation is 
small. Actually, for g„ ^ {U = 2V2), X^^iT) is con- 
stant at low temperatures as for a normal metal. Thus 
the enhancement of xC^) due to the suppression of ID 
fluctuations is less detectable in these situations. These 
numerical results suggest that reasonable predictions for 
the pressure dependence of xiT) in organic conductors 
can be made as is discussed in §5. 

4-6 Effect of nesting deviations on the susceptibility 

So far we have examined the susceptibility for the Q- 
ID systems with the perfectly- nested Fermi surface, i.e., 
t±2 — 0. Since the actual Q-ID conductors exhibits the 
imperfectly-nested Fermi surface due to the various kinds 
of interchain hopping, we examine the effect of nesting 



Fig. 17. (Color online) Magnetic susceptibility for U = 0.8t, Vi = 
0.5t, V2 = OAt, N = 15 with t^ = O.lt, tj_2 = (perfect nesting) 
and tj_2 = O.Olt (imperfect nesting). The inset shows the re- 
sponse functions as a function of the temperature. The subscripts 
i of SDW-I-CDW and dSS+/TS indicate the each response func- 
tion for tj_2 = (j = 0) and tx2 = O.Olt {i = 1). 



deviation by focusing on the simple case of a finite i j_2 . 
In Fig. 16, the magnetic susceptibility is shown as a 
function of temperature with [U,Vi,V2) = (4.0t, 0,0), 
where the solid line and the cross denote the suscep- 
tibility for the case of the perfect-nesting (t_L2 = 0) 
and that of the nesting-deviation (t_L2 = O.li^), respec- 
tively. The present case, t±2/t± =0.1, is much deviated 
from the perfect nesting since the quantity t±2 treated 
as the tight-binding approximation in eq. (5) is given 
by t±2/t± — 0.035 for quarter-filling. However, x(T) for 
t±2/t± = 0.1 and t_L2 = are essentially the same, at 
variance with the response function for SDW, Xs{T), 
shown in the inset of Fig. 16, which levels off below 
T ~ t±2, showing a clear effect of nesting deviations. 
(For the details of the suppression of the response func- 
tions and the nesting property, see Ref. 7,26) In this case 
c?-wave singlet superconductivity becomes dominant at 
low temperatures. When t±2 increases, the ground state 
moves from a SDW state to a d-wave singlet supercon- 
ducting states but x(^) remains essentially unchanged. 
Such a difference between x and Xs originates from the 
coupling of the scattering process with the momentum 
transfer between in- and out-going electrons. The con- 
tributions from electron-electron scattering at small mo- 
mentum transfer q ~ are essential for x(^)i whereas 
Xs{T) is determined by processes at large momentum 
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transfer Q ~ (2A;F,7r), for which nesting conditions have 
a strong effect. This fact can be understood easily from 
the property of the Lindhard response function Xo(q); 
where xoi/l ^kp) is very sensitive to the shape of the 
Fermi surface.'^'' 

For the case with the finite long-range interaction, we 
plot the numerical results of magnetic susceptibility and 
some response functions for U = 0.8i, Vi = 0.5f, V2 = 
QAt in Fig. 17. For such a case, V2 = U/2, and SDW co- 
exists with CDW in the ground state for t^2 = 0, whereas 
rf-wave singlet and /-wave triplet superconductivity co- 
exist for i_L2 = 0.1t_L. The difference in x(^) between 
i±2/i_L = and t±2/t± = 0.1 is extremely small in spite 
of this coexistence. 

5. Comparison with experiments 

In this section, we compare the present results for 
X{T) with those known for the Q-ID organic conduc- 
tors (TMTCF)2X (C =S, Se; X =FFe, Br). The lat- 
ter are summarized as follows.'^' ^ For the compounds 
(TMTTF)2X {X = PFe, Br), xi^) has a maximum 
at T ~250-300K, and decreases with an upward cur- 
vature as temperature decreases. On the other hand, 
for (TMTSF)2PF6, x(T) the maximum of x{T) is above 
room temperature, and the susceptibility decreases ap- 
proximately linearly with temperature. The temperature 
variation of x(T) is therefore stronger for the TMTTF 
salts. Since the present results for the Q-ID x(^) agree 
qualitatively with the ID case, one can discuss its tem- 
perature dependence in terms of the results of the ID 
system. The stronger temperature dependence of x{T) 
for the TMTTF salts suggests two possibilities: either (i) 
a large U or (ii) both a large U and Vi . The weaker tem- 
perature dependence of x{T) in the TMTSF salts would 
indicate two possibilities, namely either (i) a small U or 
(ii) a reasonable combination of U, Vi and V2 couplings. 
The charge ordered state, expected for large Vi, has been 
observed in the TMTTF salts. ^"^ Thus both U and Vi are 
expected to be sizable in this family. In the TMTSF salts, 
it is plausible that U, V\ and V2 are large, owing to the 
SDW and CDW phase coexistence and the possibility of 
triplet superconductivity for large U -\-V\-\- V2.^°'^^'^^ 

Based on the present numerical results shown in §4, we 
can suggest the following qualitative predictions about 
the pressure dependence of x(^) to make the character 
of these Q-ID conductors more clear. (Here, we assume 
that the main effect of the pressure is an increase of tj^.) 

(1) The enhancement of x{T) even at low temperature 
(T < O.lt) is essentially ascribed to small U. 

(2) The behavior of x{T) described by the enhancement 
at high temperature (T > OAt) and the reduction 
at low temperature (T < O.lt) is attributable to the 
large (7(~ At) or Vi. 

(3) The reduction of x{T) in the high temperature do- 
main (T > OAt), results from the large V2{^ U/2). 

However, from the pressure dependence of x(^)) it is still 
unclear if these conductors are characterized by interac- 
tions with large U or large ^7-1-^1. Since the behavior of 
x{T) with < t±2 < t± is almost identical to the case 
where t±2 = 0, as discussed in §4.6, the above scenarios 



deduced from t±2 = are also valid when nesting devia- 
tions are present and superconductivity is stabilized. 

6. Summary and conclusion 

In this work, we have extended the Kadanoff- Wilson 
renormalization group (RG) method to the Q-ID sys- 
tems. In addition to standard logarithmic terms of the 
one- loop RG, the present iV-chain RG scheme also in- 
cludes non logarithmic contributions of the Landau and 
finite momentum Cooper channels. The full (transverse) 
momentum dependence is retained for all contributions. 
These are essential to the description of long wavelength 
spin correlations and enter as key ingredients in the cal- 
culation of magnetic susceptibility at finite temperature. 

The ID to Q-ID crossover in the temperature (energy) 
dependence of the couplings has been obtained. When 
T > Tx ~ 0{t±), the flows of the couplings g^{{k±}) 
are essentially independent of the momenta and almost 
the same as in the ID case. When T < T^, the flows of 
9v{{k±}) are found to difi^er and deviations from the ID 
case strongly depend on the degrees of quantum inter- 
ference between various scattering channels that become 
{^±1,2,3} dependent. 

The non-perturbative influence of t± on the magnetic 
susceptibility, xi^), for Q-ID case has been calculated 
by the RG technique at the one loop-level and the re- 
sults compared to different RPA approaches. From the 
A''-chain RG, it is found that for sizable U(> At), the 
magnetic susceptibility is lower than the ID case at low 
temperature (T < O.lt), but it is enhanced by t± at high 
temperature. This contrasts with the RPA one, xkpa{T), 
which shows a t_L-enhancement in the whole temperature 
region. Such a behavior shows the significant role of t±, 
which depends on temperature. At high temperatures, 
t± suppresses the ID fluctuation, while it enhances the 
antiferromagnetic spin- fluctuation at low temperatures. 
We also note that , for small U (e.g., U < t), both RG 
and RPA always give the enhancement of x(T) by t±. 
The qualitative difference between RPA and RG essen- 
tially originates from the effect of mode-mode coupling 
terms which contribute to the RG flow. 

When long-range Coulomb repulsion is added, in par- 
ticnilar the ncxt-iiearest-neighbor repulsion V2, the tem- 
perature interval where the susceptibility decreases with 
respect to the ID limit extends to higher temperature. 
With such an effect of t±, it may become possible to es- 
timate the magnitude of long-range interactions in Q-ID 
metals through the pressure dependence of x(^)- 

We have also examined the effect of the nesting devi- 
ation on the magnetic susceptibility by introducing the 
next-nearest-neighbor transverse hopping t±2 > 0. The 
response function for CDW and SDW (the 2fcF-charge 
and spin fluctuations) saturates below T ~ t±2, while 
the eff'ect of the nesting deviation is hardly seen for the 
magnetic susceptibility and the response functions of the 
superconductivity. Namely, alteration of nesting does not 
contribute to x(^) due to the small- momentum transfer, 
in contrast to Xc.s(T) which are mainly governed by the 
large-momentum transfer (~ 2A;f) processes. 
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Appendix A: Contraction of the Landau chan- 
nel 

One of the main improvements of the present A^-chain 
RG from the previous Kadanoff- Wilson RG is that it 
can take into account the Landau channels, which arc 
non-logarithmic but give relevant contributions at finite 
temperatures. Here we show how to carry out the con- 
traction of the Landau channels in detail. 

The shell average of the quadratic term of the "Lan- 
dau" action is 



Landau contraction yields 
1 

2< 



(A-6) 

Similarly, the contraction of the Cooper -|- channel at 
the one-loop level is 
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xO:,^pJq,k,)0,,pMh)- (A-1) 

Here, the definition of Oip and 040- is the same as Op and 
Oa, respectively. The corresponding diagrams are shown 
in Fig. 4 (c). The shell average is given by 

shell 
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i.,i.'=4||,4_L{fei 3,g'}. fe' 

X A:{q',k')A,,{q',k')g^,{k'^,q'^ - k'^MrA)) 
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The corresponding diagrams are shown in Fig. 4 (d). The 
shell average is 

shell 
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Ml, Q^, fcl) = vpq + a(fcl) - U{k'± - q±), (A-4) 

vpq = vpiks - ki) 

= U{ki±)-Uiki±+q±). (A-5) 

The function /l indicates a deviation from the logarith- 
mic divergence, but it does not reach the unity, i.e., the 
Landau channel is always less logarithmic. Finally, the 



Ac+{q',q'^,k'^) = Vpq' + U{k±) + U{k± - q±), 

(AlO) 

vpq' = vpiki + k2 - 2kF) 

= -U{ki±)-i±{q'±-ki±). (A-ll) 
The Cooper-h contraction yields 
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The contributions of the Cooper and Peierls channels 
bubbles /c,p and that of Landau and Cooper+ channels, 
Jl,c+! are approximated for small t± as 
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Fig. B-1. Temperature dependences of the functions /(T) (solid 
line) and g{T) (dashed line). 
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Then, the flow equations (52) can be rewritten in the 
form: 
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From these equations, the couplings in Q-ID on the 
2T J basis of a perturbation expansion in t± are obtained. 

The couplings for q = are shown in §4. It is use- 
ful to define the functions f{T) = j'dl f{Eo{t,/AT), 
g(T) = Jd£ g{Eo{£)/4T). The temperature dependences 
of f{T) and g{T) are shown in Fig. B-1. In the limit 
of T ^ 0, each function is calculated analytically as 
Aciqc,qc±,k'^)]^ f{T) 7C(3)/47r2 ~ 0.213, g{T) 0. As is seen from 
J Fig. B-1, /(T) and g{T) can be treated as constants at 
low temperatures (T < 0.2t). 

•y^ I Ap{Q, Q±, kjj I Appendix C: Response functions for staggered 
1./ I ^-^ J density-wave and superconductiv- 

ity 



V fr r , 9{Eoie)/^T) ^ 

-Uo-(U4p - 2^ 

fc', 



2T 




Ap{Q,Qi_ 




2T 




f Ai,{q,q± 




\ 2T 






(B-5) 



C.l SDW and CD W responses 
We can calculate the SDW and CDW response func- 



s';,[r,^]^J2^h,iQ)o;iQ,k) + c.c:_ 



(C-l) 



where ^ = c, s. Here, we concentrate on the response to 
the source field with the nesting vector Qo = {2kp^n). 
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This commensurate nesting vector is appropriate for 
the present Fermi surface, ahhough the incommensurate 
nesting vector is needed for an arbitrary Fermi surface. 
The one-loop corrections yield 

+ XMW/iM(Qo)/iM(Oo), (C-2) 
The pair vertex part is expanded in the form 
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Fig. C-1. Response functions of the SDW (solid lines) and 
the d-wave singlet superconductivity (dashed lines) for t_i_ = 
^ 0,0.01,0.1* and U = 3t. 



In the present paper, we employ the approximation 
gti{k±,Qo± - fcj., ^ gfiik±,Qo± -k'^,k'j_) and carry 
out the summation over k'j_ first, since we consider the 
full-gapped density wave state so that the fe'^-dependence 
is irrelevant. The response function in units of {TrvF)~^ 
is given by 
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C.2 Superconductivity response 

The quasi-one-dimensionality, i.e., the use of a finite- 
t± , makes anisotropic Cooper pairing possible such as d- 
wave singlet or /-wave triplet, which cannot be realized 
in pure ID systems. Thus we consider not only the full- 
gapped superconductivity, which have been covered by 
previous RG calculations, but also for anisotropic super- 
conductivity corresponding to an order parameter with 
nodes. We add the 5^' to the action: 

S'^^'ir,^] = 5][V(9c)A;,(gc,^) +C.C.], (C-6) 
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where /i' — ss, ts. The one-loop corrections yield 
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At the one-loop level, the pair vertex part is expanded 
in the form 

zl::\£+di)^z^:\£) 

-z^;\e)d;hc{qco,qco±,0)di. (C-8) 
The flow equation for zf?-* is 
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^l^?-* ) is the Fourier coefficients of the 



where C;,/(= , 

coupling gss,ts, which are defined in the form 
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(For the details of the pairing symmetry, see Ref. 21.) 
The response function in imits of {TrvF)~^ is given by 
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The response functions of SDW and dSS for t±_ = 
0,0.01, 0.lt,U = 3t are shown in Fig. C-1. For the ID 
case, Xs{T) ^ T^°' [a ~ Gp/2), which fully agrees with 
the well know results of the previous RG.^^ For ti_ > 0, 
XsiT) is enhanced below Tx. The response function for 
dSS, Xss^ does not grow for ID {tj_ — 0), because the 
pairing with nodes is very weak due to the geometric re- 
striction of ID or due to lack of transverse single coher- 
ence. On the other hand, in Q-ID systems, such a pairing 
becomes possible due to the warping of the Fermi sur- 
face. As a result, Xss\t) rapidly grows below Tx, though 
it is less singular than Xs{T) in the case of a perfectly 
nested Fermi surface. 
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